#include <gmlib/sm9/internal/sm9_field.h>
#include <stdexcept>
#include <cstring>

using namespace sm9::internal;

void test_sm9_fp_sqrt()
{
    sm9_fp_t    a, c;

    static std::uint8_t a0[32] = {0x0a,0xcf,0xf0,0xd7,0x69,0x5f,0x43,0x85,0x57,0x61,0x30,0xde,0xc8,0x7e,0x0a,0x75,0x34,0x16,0xa2,0xb1,0x7a,0x7e,0xa5,0x2d,0xd5,0x0f,0xf9,0x48,0xbe,0x85,0x4d,0x1a};
    static std::uint8_t c0[32] = {0x30,0x8b,0x2b,0x22,0xf3,0x7a,0xf1,0x14,0xd7,0xf1,0xb7,0x64,0xd8,0xef,0x87,0x25,0xfe,0xcd,0xec,0x21,0x90,0x74,0xab,0x5a,0x25,0xd4,0x6c,0x3d,0x85,0x56,0x7c,0x7c};
    sm9_fp_from_bytes(a, a0);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a1[32] = {0x0e,0x93,0x0c,0x32,0x3f,0xc9,0x00,0x09,0x8a,0x0d,0xea,0xe3,0xc8,0xae,0xf9,0x73,0x63,0xb1,0x23,0xae,0x08,0x4f,0x5e,0x52,0xf7,0xb5,0x69,0x17,0x54,0x99,0xd7,0x05};
    static std::uint8_t c1[32] = {0x76,0x3b,0x2b,0x52,0x7d,0x36,0x2c,0xf0,0x99,0x1e,0x05,0xbb,0xb4,0xc4,0x8f,0x23,0x0b,0x5a,0x2d,0xb0,0xd9,0x1f,0xa2,0xe1,0x73,0xbe,0x8c,0x16,0xd4,0x19,0x6c,0x85};
    sm9_fp_from_bytes(a, a1);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a2[32] = {0x81,0x98,0x27,0xa2,0xf2,0x0c,0x75,0x22,0x1e,0x25,0x84,0x30,0xe4,0x30,0xad,0x1d,0x68,0x49,0xa9,0x8c,0x4f,0x95,0xef,0xaa,0x52,0xdf,0x33,0x61,0x56,0x4a,0x26,0x55};
    static std::uint8_t c2[32] = {0x42,0x2c,0xc5,0xd3,0x28,0xc9,0x8b,0xa2,0x56,0x92,0x0b,0xe1,0xca,0x41,0xd3,0x7f,0xc8,0x1d,0xfa,0x55,0x03,0x4d,0x70,0x39,0x62,0x31,0x92,0xec,0x66,0x72,0xe6,0x87};
    sm9_fp_from_bytes(a, a2);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a3[32] = {0x33,0xd2,0x12,0x04,0x04,0xe4,0x5d,0x30,0x62,0x04,0x98,0x19,0x02,0x65,0x40,0x10,0x4b,0x64,0x3a,0xfe,0x0d,0x53,0xfa,0x48,0xa4,0xf5,0x6a,0xb3,0xc7,0x97,0x6c,0x05};
    static std::uint8_t c3[32] = {0x9c,0xfb,0xb0,0xcf,0xb2,0x41,0xb5,0x56,0xd9,0x99,0xd3,0x71,0xd3,0x6c,0xb3,0x02,0xd2,0xd0,0x4d,0x49,0xc9,0x7a,0x09,0xb8,0x2a,0x22,0x33,0x89,0x2e,0x4b,0x22,0x35};
    sm9_fp_from_bytes(a, a3);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a4[32] = {0x75,0x9d,0x77,0x69,0x0f,0xed,0xb1,0x76,0x6d,0xd0,0xe7,0x22,0x68,0xfa,0x34,0x4d,0x63,0xc5,0x76,0xf4,0x91,0x75,0x3b,0xda,0x4c,0x92,0xe4,0x8c,0x05,0x48,0x58,0x9f};
    static std::uint8_t c4[32] = {0x3b,0x1b,0x34,0xf8,0x12,0xda,0x58,0x9f,0xd4,0xb6,0xd5,0xff,0x4f,0x10,0x76,0x06,0x4a,0x5d,0xab,0x56,0xa8,0x02,0x51,0xbf,0xad,0xcc,0x2a,0x5d,0xe2,0x4a,0x67,0x8b};
    sm9_fp_from_bytes(a, a4);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a5[32] = {0xa3,0xe7,0xe2,0x0b,0xc3,0xf9,0x54,0x1a,0x88,0x81,0xc8,0xeb,0xe9,0x3d,0x1f,0xc5,0xad,0x2d,0xae,0x37,0x58,0xa3,0xab,0x42,0xc7,0xee,0x4c,0xf1,0x92,0xcc,0x69,0x8f};
    static std::uint8_t c5[32] = {0x6d,0xc3,0xa8,0xc0,0x53,0x82,0x28,0x4b,0xd1,0x68,0xa4,0xef,0x58,0xf8,0x34,0x5c,0x0a,0x70,0x4e,0xe4,0xf2,0x54,0xf4,0x70,0x4c,0xcf,0x6d,0x1b,0xb6,0x20,0x47,0xeb};
    sm9_fp_from_bytes(a, a5);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a6[32] = {0xa8,0x1f,0xae,0x1b,0x35,0xfc,0x9b,0x6e,0x0b,0xc9,0xfd,0x69,0x36,0xea,0x65,0x8d,0xe9,0x2c,0x58,0xb5,0x15,0x79,0x04,0xe3,0x5d,0x52,0xce,0x00,0x75,0x99,0x1c,0x1c};
    static std::uint8_t c6[32] = {0x51,0xc5,0x7d,0xc6,0x46,0x79,0x1c,0x90,0x43,0x6f,0x51,0x76,0x92,0x02,0x77,0x16,0xb6,0x2f,0xd3,0x0e,0xb0,0xba,0xe2,0xac,0x5f,0xaa,0x84,0x9f,0x14,0xf2,0xc9,0x80};
    sm9_fp_from_bytes(a, a6);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a7[32] = {0x20,0xd6,0xa3,0x3b,0xf8,0xe3,0xb7,0x2e,0x44,0xf3,0x4f,0xb3,0x85,0x8b,0x0a,0x1e,0x88,0x8f,0xee,0x9e,0x93,0xa9,0x10,0xd9,0xc1,0xdc,0x38,0xc3,0xc3,0x48,0x0d,0x14};
    static std::uint8_t c7[32] = {0x3a,0xd2,0x5a,0x05,0xfd,0x39,0xa2,0x07,0x36,0xe1,0xdc,0xf4,0x9e,0x61,0x81,0x4c,0x65,0x4d,0x3f,0xcd,0xbc,0x54,0xe3,0xae,0xbf,0xac,0xb0,0xdf,0x6a,0xb6,0x0f,0xef};
    sm9_fp_from_bytes(a, a7);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a8[32] = {0x5e,0x09,0x02,0xcb,0xa7,0x3b,0x4b,0x8d,0x9b,0xbe,0xd4,0xc2,0x9d,0x9e,0xa9,0x58,0x5e,0x31,0x49,0x2c,0x4b,0x71,0x6b,0x77,0x7a,0xec,0x60,0x92,0xe2,0xfe,0x0e,0xd4};
    static std::uint8_t c8[32] = {0xa7,0xde,0x64,0x84,0x8e,0xdf,0x61,0xb4,0x81,0x8d,0x0b,0xa1,0x04,0x7a,0x2c,0x26,0xc3,0xb4,0x81,0xff,0x3d,0x45,0x81,0xb7,0xc7,0x91,0x00,0x3a,0xcd,0xcb,0x14,0x06};
    sm9_fp_from_bytes(a, a8);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a9[32] = {0x8c,0x73,0x4d,0x1e,0xc3,0x2a,0xe3,0x63,0xfb,0xf0,0xd1,0xd5,0xc5,0xca,0x91,0x73,0xf0,0x80,0x3e,0xbf,0x7c,0xc7,0xd5,0xee,0x6b,0x8e,0xd4,0x3f,0xe0,0x92,0x59,0x04};
    static std::uint8_t c9[32] = {0x9a,0xb3,0x58,0xcc,0x38,0xd2,0x3e,0x63,0x4b,0xe1,0x80,0x4e,0x87,0xa3,0x04,0xbb,0x96,0x37,0x91,0x2e,0xab,0xef,0x23,0xf1,0xf9,0xc2,0x10,0x1c,0xba,0xb7,0x9f,0x2e};
    sm9_fp_from_bytes(a, a9);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a10[32] = {0x0b,0xd0,0x48,0x19,0x3e,0x53,0x18,0x85,0xa0,0xfd,0xc5,0x7d,0x2b,0x19,0x63,0x62,0x15,0x21,0x8d,0x3d,0x50,0x97,0xda,0x4b,0x2b,0x97,0x8a,0x71,0xba,0x53,0xa6,0x80};
    static std::uint8_t c10[32] = {0x86,0x10,0xd6,0x9c,0xdf,0x8a,0x30,0x49,0xd0,0xf5,0xe1,0x83,0x05,0x5e,0x98,0x3e,0x43,0x78,0xc7,0x23,0xfe,0x3f,0xc8,0x65,0x49,0xea,0x94,0xf3,0x70,0x9f,0x5b,0x29};
    sm9_fp_from_bytes(a, a10);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a11[32] = {0x87,0xc7,0xda,0x00,0x5c,0x8b,0xe1,0x3e,0x02,0x25,0xe1,0x2b,0x7a,0xec,0x9f,0xb9,0x1d,0x25,0x69,0x68,0x7e,0xe8,0x57,0xad,0xd0,0xc7,0x36,0x3a,0x87,0x83,0x77,0xee};
    static std::uint8_t c11[32] = {0xb0,0xb9,0xa4,0x35,0x0e,0xb2,0xe1,0xb4,0x8e,0x66,0x42,0x99,0xe5,0x39,0xa8,0xec,0x8d,0x8d,0x64,0x8d,0x3a,0xec,0x8b,0x01,0xaa,0x09,0x31,0xa6,0x29,0x27,0x35,0xcf};
    sm9_fp_from_bytes(a, a11);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a12[32] = {0x39,0x28,0x35,0x66,0x82,0xbd,0xeb,0xf8,0x27,0xe2,0x47,0xca,0x7c,0x28,0xc7,0x3d,0x21,0x46,0x51,0x9f,0x22,0x73,0x7f,0x3b,0xf5,0x2c,0x9c,0x40,0xf7,0x4e,0x14,0x81};
    static std::uint8_t c12[32] = {0x57,0x84,0x93,0xeb,0xcc,0x7f,0xc4,0x5d,0x1a,0x46,0x5c,0xa6,0xcb,0xae,0x05,0x8a,0x81,0xce,0xe0,0xf5,0x71,0xb2,0x70,0x14,0x21,0xc5,0xdf,0x70,0xac,0xb7,0xa2,0x04};
    sm9_fp_from_bytes(a, a12);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a13[32] = {0xa4,0x35,0xbd,0xba,0x86,0xaa,0x8a,0xcb,0x4a,0x50,0xcb,0xfb,0xd7,0x8f,0xe2,0xa9,0x9c,0xcc,0x6a,0x65,0x29,0x1a,0xf8,0x30,0x84,0xea,0x25,0xc7,0xb9,0x96,0xa2,0x94};
    static std::uint8_t c13[32] = {0x90,0xca,0x33,0x29,0xd8,0xe3,0x9c,0xfd,0x63,0xc1,0xaa,0x98,0x13,0x59,0xc0,0x92,0xf7,0xc2,0xbd,0x66,0x0a,0xee,0x87,0xff,0xc0,0x7a,0x86,0x10,0x4a,0xd6,0xa8,0x91};
    sm9_fp_from_bytes(a, a13);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a14[32] = {0x18,0x6f,0x4e,0x30,0x78,0x6c,0x3d,0xf0,0x47,0x83,0xcf,0xae,0xb8,0xc7,0x7f,0xd2,0xe0,0x86,0xf0,0xf6,0x8c,0x92,0x9e,0xc5,0x33,0x46,0x59,0x11,0x15,0x43,0x9b,0x53};
    static std::uint8_t c14[32] = {0x52,0x3a,0x28,0xf1,0xac,0x4a,0x58,0x71,0xa0,0xae,0x68,0x71,0x44,0xce,0x5c,0xba,0xa4,0x82,0xf5,0x34,0xf5,0x48,0x46,0xe9,0xf0,0x16,0x8c,0x0a,0x22,0x6f,0x94,0x77};
    sm9_fp_from_bytes(a, a14);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a15[32] = {0x05,0xdd,0x89,0x2e,0x53,0x0c,0x52,0x28,0x60,0xd3,0x71,0x2e,0x0e,0x07,0x63,0x6f,0x9f,0xb0,0x89,0xdb,0x2a,0x08,0xde,0xe4,0x5b,0x58,0x01,0x40,0x4e,0x5f,0x1c,0x41};
    static std::uint8_t c15[32] = {0x97,0xdb,0xf6,0x57,0x94,0x72,0xe3,0x71,0x3e,0x72,0x75,0xd8,0x0a,0xf4,0x80,0xbe,0x70,0xdd,0x2a,0x05,0x36,0x17,0x3c,0x8c,0xf4,0x86,0x6d,0x59,0xd5,0xba,0xd2,0xaf};
    sm9_fp_from_bytes(a, a15);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a16[32] = {0x64,0x01,0xaa,0x21,0x26,0xbc,0xa8,0x04,0xff,0xed,0x02,0x85,0xe3,0x15,0xb1,0xe4,0xa4,0x1b,0x48,0xd4,0x78,0x93,0x5f,0x97,0x69,0xfe,0xa2,0x44,0x35,0x7e,0x37,0x29};
    static std::uint8_t c16[32] = {0x29,0xaf,0x82,0xdc,0x82,0xa4,0xcf,0xad,0xa8,0xe1,0xce,0xbb,0xfe,0xa2,0x21,0xf2,0xdc,0x12,0xe5,0x8f,0x33,0xd9,0xcb,0x53,0x4a,0x81,0x49,0x51,0x38,0x68,0x0c,0x40};
    sm9_fp_from_bytes(a, a16);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a17[32] = {0xa8,0x35,0x3b,0xde,0x65,0xdf,0x8a,0xad,0xc5,0x73,0xee,0xd5,0xc3,0xec,0x82,0x4b,0x35,0x51,0x10,0xb0,0x14,0x0e,0x6a,0xed,0xf4,0x0f,0xa1,0x03,0x1f,0xd9,0x43,0x81};
    static std::uint8_t c17[32] = {0x14,0x0b,0x9c,0xf1,0xab,0xe7,0x0e,0xed,0x82,0x0a,0x90,0x47,0xfb,0x99,0xf7,0x4d,0x4f,0x78,0xc3,0x59,0x17,0xc8,0xa5,0xaa,0x74,0x2d,0xa6,0x20,0xa5,0x25,0x8c,0x47};
    sm9_fp_from_bytes(a, a17);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a18[32] = {0x6f,0x54,0x67,0x23,0xa8,0xb4,0xeb,0x75,0x34,0x72,0xa3,0x0b,0xb2,0xcd,0x53,0x0c,0xdb,0x26,0xb5,0x43,0xf0,0x9d,0x76,0x62,0x55,0x6b,0x87,0xa3,0xc1,0xcb,0x8a,0xc3};
    static std::uint8_t c18[32] = {0x65,0xf5,0x2e,0xfa,0xdf,0xdd,0x1c,0x0e,0xa4,0x45,0x78,0xac,0x52,0x8b,0x60,0x31,0xc6,0x17,0x4c,0x5a,0x19,0xcb,0x22,0x6f,0xae,0xce,0xb2,0x3b,0xc1,0xd4,0x08,0x7d};
    sm9_fp_from_bytes(a, a18);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a19[32] = {0x48,0x3e,0xd8,0x19,0x3a,0x3f,0x32,0xbb,0x8f,0xfe,0x1b,0x9c,0xd1,0x57,0xa1,0x0e,0xe6,0x8d,0xde,0x43,0xb6,0x5e,0x3f,0xf4,0x44,0x63,0x99,0xe5,0x1e,0x77,0xbc,0x0e};
    static std::uint8_t c19[32] = {0x09,0xd7,0xf1,0x26,0x04,0x39,0x41,0x98,0xe8,0xbb,0x0f,0x68,0x4a,0xa0,0x1f,0xff,0x30,0x3b,0xd9,0xda,0x0f,0xd8,0xa2,0x57,0xda,0x0b,0x7a,0xab,0x85,0xfa,0xfc,0xda};
    sm9_fp_from_bytes(a, a19);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a20[32] = {0x35,0xf6,0x06,0xcd,0x8b,0xb2,0x15,0x48,0x98,0x6d,0x22,0xf5,0xd5,0x7f,0x85,0x7b,0x72,0x7d,0x79,0xb6,0x90,0x04,0xa6,0xed,0x98,0x79,0xbc,0x19,0x94,0xe2,0x4e,0x62};
    static std::uint8_t c20[32] = {0x9b,0xca,0xeb,0x26,0xa0,0x78,0xe4,0x7a,0xcd,0x12,0x6d,0x92,0x54,0xcf,0x6e,0x3d,0xf8,0xef,0xcc,0xab,0x00,0xae,0xc9,0x95,0xbd,0xdf,0xa0,0x13,0x25,0x4f,0xf9,0xc7};
    sm9_fp_from_bytes(a, a20);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a21[32] = {0x5e,0xa7,0xa3,0x05,0x83,0x6b,0x88,0xd3,0x87,0x15,0xb2,0x31,0x7d,0x8d,0x11,0xa0,0xd7,0x95,0x48,0xa3,0xd2,0x7b,0x69,0xd6,0x52,0x02,0x83,0x65,0xbb,0x17,0x41,0x52};
    static std::uint8_t c21[32] = {0x62,0x09,0xa2,0xb4,0xbe,0x8d,0x72,0x80,0xa5,0x21,0x02,0x9d,0x38,0xa7,0x3e,0x9e,0x97,0x15,0xc4,0x03,0xef,0x4c,0x48,0x5f,0xa9,0x85,0xf6,0x95,0xdf,0x8a,0x72,0x1d};
    sm9_fp_from_bytes(a, a21);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a22[32] = {0x2b,0x4f,0x11,0xe2,0x00,0x64,0xf6,0x91,0x15,0x3f,0x3c,0x4b,0xd6,0xe5,0xc0,0x21,0x2f,0x17,0xe0,0x7e,0xf3,0x74,0x96,0xa1,0xfe,0x60,0xda,0x64,0x55,0x7f,0x25,0x4f};
    static std::uint8_t c22[32] = {0x26,0x78,0x7c,0xe8,0x11,0x0b,0x89,0xa9,0x42,0x82,0x59,0xa5,0x26,0xb4,0xc8,0xb4,0xdd,0x18,0xe7,0xbb,0x3e,0x8c,0xa4,0x41,0x1f,0xd3,0x9b,0xe8,0x00,0x7c,0xf9,0x4a};
    sm9_fp_from_bytes(a, a22);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a23[32] = {0x71,0xf2,0x91,0x3e,0x38,0xe1,0x0d,0x57,0xa2,0x50,0x38,0x63,0xcb,0xa4,0x85,0x41,0x6e,0x27,0x58,0x99,0x25,0x65,0x51,0xfb,0x4d,0x45,0xe8,0xaf,0x75,0x72,0x70,0x38};
    static std::uint8_t c23[32] = {0x1d,0xa6,0xfd,0xda,0xfe,0x0c,0xfd,0xa4,0x80,0x83,0xb3,0xdb,0x5f,0xf6,0xa9,0x32,0x82,0xd0,0x6c,0x62,0x6e,0x1f,0x34,0x73,0xbc,0x02,0x8e,0xa3,0xef,0xb8,0xc4,0x6e};
    sm9_fp_from_bytes(a, a23);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a24[32] = {0x25,0x4d,0x88,0x1f,0x9e,0x7a,0x52,0x88,0x9d,0x42,0x4c,0xf0,0xe4,0x55,0xcc,0x09,0xf3,0xb2,0x75,0x29,0x0a,0xf7,0xa9,0x94,0xa3,0xfe,0x52,0x7f,0x2f,0xaa,0x76,0x07};
    static std::uint8_t c24[32] = {0x24,0xe1,0x98,0xf5,0x00,0x28,0xd4,0xa1,0xc1,0x90,0x50,0xea,0x70,0x71,0xb3,0x3b,0xbc,0x93,0x23,0x9d,0x10,0x77,0xef,0x5d,0x2f,0x6f,0x0f,0x59,0x4a,0xb3,0xa0,0xc6};
    sm9_fp_from_bytes(a, a24);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a25[32] = {0x15,0x2f,0x0d,0x59,0x5e,0xf7,0xa4,0x3d,0x78,0xbc,0xbf,0x9f,0x92,0x52,0xa9,0xd4,0xb3,0x17,0x51,0x44,0xc7,0x8f,0x46,0x53,0x92,0xbd,0x49,0xf6,0x85,0x02,0x29,0xf2};
    static std::uint8_t c25[32] = {0xa9,0xb3,0x25,0x27,0x78,0x93,0x91,0x1d,0xb7,0xb5,0xd7,0x77,0x61,0x14,0xbc,0x0a,0x5c,0x64,0x18,0xdf,0x59,0x4c,0x23,0x87,0xcd,0xf0,0xac,0x99,0xa3,0x6c,0x9a,0xba};
    sm9_fp_from_bytes(a, a25);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a26[32] = {0x95,0x11,0xf3,0x11,0x78,0x0c,0xb7,0xdd,0x73,0xd0,0xc0,0xee,0x49,0x82,0x84,0xb0,0xb6,0x8b,0x03,0xa2,0x86,0xfe,0x0e,0x4e,0x1a,0x9f,0xd3,0x9c,0xbd,0x69,0x80,0x22};
    static std::uint8_t c26[32] = {0x54,0x5e,0x44,0x98,0xa2,0x85,0x1f,0x75,0x6b,0xdb,0x33,0x92,0x7a,0x0c,0x9d,0x97,0xc9,0x28,0x61,0x1a,0x29,0xdf,0xa9,0x3b,0xef,0x22,0xe0,0xc1,0x38,0x54,0xde,0xf2};
    sm9_fp_from_bytes(a, a26);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a27[32] = {0x3c,0x00,0xa2,0xfc,0xfa,0x44,0x66,0xb1,0x51,0xee,0x9a,0xe4,0xad,0xca,0xba,0xa5,0x50,0xfd,0xba,0xcf,0x84,0x1a,0x35,0x52,0x0f,0x9e,0x23,0x29,0x10,0x42,0x59,0x9c};
    static std::uint8_t c27[32] = {0x0c,0x65,0xf0,0x8f,0x66,0x94,0xdc,0x75,0x3c,0x15,0x8c,0x9a,0xb4,0x6e,0xc7,0x44,0x73,0x44,0xc8,0x8e,0x79,0x89,0x9f,0xac,0x14,0x57,0x2d,0x38,0x16,0x2c,0x5e,0x61};
    sm9_fp_from_bytes(a, a27);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a28[32] = {0x5e,0xf7,0x64,0x45,0xc9,0x69,0x08,0xa9,0xad,0x7e,0x73,0x26,0x68,0xaa,0xeb,0x8b,0x34,0x88,0x6a,0xa5,0x66,0x09,0x51,0x12,0xdd,0xac,0x1a,0x52,0xaa,0xb8,0xef,0xe4};
    static std::uint8_t c28[32] = {0x07,0x1d,0xaf,0x58,0x76,0x4a,0x58,0xfc,0xee,0x78,0xa3,0x95,0xc9,0xfb,0xa3,0xd6,0x24,0x8e,0x57,0x47,0x09,0xf7,0x2f,0x8f,0x09,0x43,0xac,0xe0,0x68,0x50,0xa2,0xea};
    sm9_fp_from_bytes(a, a28);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a29[32] = {0x29,0xb9,0xc1,0x7c,0x9a,0x8b,0x88,0x6d,0x40,0x35,0x5f,0xea,0xc4,0x14,0x80,0x39,0xa3,0xba,0x2f,0x1c,0x3a,0x7e,0x50,0x90,0x12,0x1b,0x94,0x9d,0x68,0x03,0xce,0x66};
    static std::uint8_t c29[32] = {0x2d,0xf9,0x8f,0x11,0xde,0xe5,0xa3,0x06,0x8c,0x9b,0x69,0x15,0x2b,0x45,0xc3,0x41,0xc1,0xd6,0xee,0xd5,0x4f,0x3a,0x83,0x32,0x02,0xb3,0x9f,0xbb,0xe1,0x49,0x62,0xc5};
    sm9_fp_from_bytes(a, a29);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a30[32] = {0x5b,0xf0,0x3e,0x76,0x9c,0xc1,0x8f,0x8d,0xe7,0x0b,0x30,0xbd,0xe9,0x7b,0x01,0x06,0xef,0x78,0x22,0xaa,0xfb,0xad,0xe6,0x18,0x42,0x78,0xef,0xf4,0x6f,0x72,0xd3,0x58};
    static std::uint8_t c30[32] = {0x60,0xa2,0xb4,0x1e,0x3b,0xab,0xd8,0xd4,0x71,0xa3,0x55,0xb4,0xfe,0xae,0xc8,0xb2,0xe8,0x63,0x1c,0x87,0x94,0xd8,0x13,0xd8,0x96,0x6d,0xa6,0x19,0x5d,0x49,0x9c,0xfb};
    sm9_fp_from_bytes(a, a30);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a31[32] = {0x1c,0x2f,0x92,0xb0,0x77,0xc4,0xb6,0x19,0x61,0xae,0xc5,0x70,0x2e,0xb5,0x55,0xf6,0x88,0xb3,0x61,0x58,0x20,0xb2,0xff,0xd3,0x74,0x7a,0x31,0xe8,0xf3,0xa0,0x57,0x18};
    static std::uint8_t c31[32] = {0x67,0x91,0xd5,0x8b,0xd5,0x19,0x45,0xd2,0x70,0xc5,0x69,0xde,0xbc,0x34,0x94,0xd4,0x30,0x44,0x52,0x14,0x83,0xf3,0x08,0x0a,0x2a,0x0a,0xd5,0x10,0xc2,0xe9,0x63,0xe7};
    sm9_fp_from_bytes(a, a31);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a32[32] = {0xa6,0xf6,0x41,0x53,0x51,0x16,0xaf,0xee,0x02,0xf8,0x0a,0xdd,0xc1,0x4e,0xbb,0xd6,0x5c,0x9e,0x90,0x3b,0xb1,0xbd,0x74,0xe6,0x5e,0x0e,0x6c,0x5c,0x5c,0x54,0xe4,0x7c};
    static std::uint8_t c32[32] = {0xb0,0xa4,0x63,0x3d,0xf0,0x92,0x24,0xb3,0x05,0x88,0x35,0x93,0xde,0xd4,0x45,0xb1,0x7c,0xb6,0x3e,0xa2,0x53,0xb8,0x53,0xff,0x3b,0x0d,0x2c,0xa5,0xba,0x65,0xa6,0x01};
    sm9_fp_from_bytes(a, a32);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a33[32] = {0x18,0xe7,0x9c,0xc0,0xfc,0x2f,0x48,0x29,0x70,0xc1,0xc4,0x7e,0xa3,0x1e,0x29,0x8a,0xdc,0x4e,0x9c,0x25,0x8d,0xd6,0xfe,0x74,0x85,0x08,0x97,0x93,0x21,0x1b,0x9d,0x04};
    static std::uint8_t c33[32] = {0x87,0x2c,0xce,0x69,0xc6,0xb6,0x6b,0x29,0x4a,0xbb,0x55,0x21,0x73,0x62,0x48,0xba,0xd7,0x66,0xe7,0xf6,0xc8,0xa3,0x36,0x0d,0x24,0xc8,0x66,0x00,0x6d,0xd4,0x53,0xd0};
    sm9_fp_from_bytes(a, a33);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a34[32] = {0x4f,0x72,0x6e,0x66,0xd5,0x43,0xde,0x7c,0xae,0x2a,0xe2,0x10,0xef,0x9a,0x73,0x3c,0xf2,0x14,0xb2,0xcb,0xe4,0xad,0x48,0xc8,0xc1,0x18,0xb5,0xaa,0xa1,0x5f,0x62,0xa8};
    static std::uint8_t c34[32] = {0x27,0x38,0x31,0x2b,0x22,0x38,0x52,0x92,0x97,0x14,0xf9,0x78,0x26,0x89,0xa2,0x1c,0x9b,0xd0,0x17,0x62,0xa0,0xfc,0x95,0x64,0xf5,0xa0,0x1a,0x44,0x89,0x78,0xed,0xcf};
    sm9_fp_from_bytes(a, a34);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a35[32] = {0x8e,0xe6,0x80,0x6c,0xe1,0xff,0xea,0x90,0xa5,0xa4,0x15,0x78,0x3c,0x88,0xc6,0x1b,0x03,0x15,0x27,0xe3,0x52,0xcf,0x97,0x11,0x83,0xe6,0xa2,0x8f,0x4f,0xdc,0x8e,0x56};
    static std::uint8_t c35[32] = {0x9b,0x93,0xfb,0xd6,0x7b,0xb9,0x33,0xe3,0xee,0xb0,0x73,0x12,0x01,0x05,0xce,0x29,0x62,0xe1,0xda,0x43,0xf9,0x33,0xbf,0x16,0x32,0x77,0x26,0xef,0x8d,0x17,0xd3,0xec};
    sm9_fp_from_bytes(a, a35);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a36[32] = {0xb5,0xd8,0x81,0x44,0x80,0xee,0xf5,0x07,0xe5,0x76,0xed,0x8d,0x86,0x42,0x17,0xc9,0x2d,0x90,0x28,0x59,0xe4,0x1e,0x69,0x45,0xdd,0xbe,0x99,0x7e,0x19,0xb3,0xab,0xae};
    static std::uint8_t c36[32] = {0x14,0xde,0x2e,0x0a,0x82,0x6b,0xa0,0xb6,0x06,0x8f,0xa6,0x7e,0xbc,0x19,0xa5,0xe1,0xc8,0x6d,0xe6,0x27,0x9f,0x39,0x90,0xae,0x38,0xa9,0x04,0x07,0x1c,0x81,0x65,0xe4};
    sm9_fp_from_bytes(a, a36);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a37[32] = {0x95,0xe3,0x23,0x1d,0xda,0x60,0xf8,0x49,0xdf,0xf4,0x46,0x96,0x7b,0x2e,0x72,0xbf,0xc5,0xb6,0x9f,0x3f,0x65,0x75,0xa9,0x8b,0x96,0x7e,0x86,0xcb,0xdb,0x54,0xa9,0x68};
    static std::uint8_t c37[32] = {0x70,0x04,0x90,0x5e,0xa1,0x5b,0xfd,0x31,0xae,0x3f,0x9b,0x95,0x15,0x2f,0x5b,0x8b,0xc8,0x7b,0xef,0x51,0x59,0xa2,0xe0,0xa7,0x69,0xcf,0xe9,0xc9,0xcb,0x63,0x67,0xe0};
    sm9_fp_from_bytes(a, a37);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a38[32] = {0x85,0xfe,0xb4,0xbd,0x05,0x19,0x5a,0xd1,0x1b,0xe5,0x1a,0xcb,0x9b,0x09,0xb5,0x12,0xac,0x67,0x5c,0x48,0x80,0x68,0xde,0xbb,0x93,0x73,0x15,0xbb,0x09,0x64,0xfa,0x69};
    static std::uint8_t c38[32] = {0x3b,0xf0,0xec,0x2d,0x52,0x97,0xca,0xe2,0x19,0x76,0x65,0x8f,0x82,0xb1,0x7e,0xd2,0xa4,0xca,0xad,0x3a,0x90,0x46,0xf8,0x3c,0xe9,0x45,0x72,0x6f,0x05,0x21,0xe1,0x8f};
    sm9_fp_from_bytes(a, a38);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a39[32] = {0x50,0x98,0x67,0xf6,0x66,0xce,0x9c,0x1a,0xf3,0x9b,0xf3,0x0a,0xe9,0x9a,0xb8,0x97,0x62,0x50,0x65,0x8c,0xa6,0x05,0x33,0xea,0x0d,0x3f,0xfb,0xca,0x56,0x53,0xe3,0xfb};
    static std::uint8_t c39[32] = {0x20,0xd9,0x89,0x1b,0x0b,0x77,0x3c,0x4b,0x50,0x93,0x2a,0xf6,0xac,0xb7,0xdd,0xf2,0x05,0x40,0x64,0x0b,0x9a,0x12,0x7b,0xf4,0x17,0x2a,0x3b,0x28,0x9d,0x54,0x36,0xad};
    sm9_fp_from_bytes(a, a39);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a40[32] = {0x39,0xa8,0xae,0xcb,0x8a,0x60,0x06,0x4f,0x7a,0x7f,0x98,0xdb,0xff,0xdf,0x6f,0xe9,0xdd,0x82,0xf7,0xe7,0x7e,0xcf,0x77,0xca,0xa5,0xff,0xf8,0xa1,0x39,0x34,0xcf,0x57};
    static std::uint8_t c40[32] = {0x5a,0xf5,0x84,0x3b,0x06,0xc9,0x4a,0x6b,0xae,0x68,0x5b,0x25,0x5d,0x8b,0xc6,0x31,0x42,0x4e,0x6e,0x09,0x11,0x44,0x16,0xaf,0xde,0xc5,0x01,0x71,0x3d,0x64,0x76,0x42};
    sm9_fp_from_bytes(a, a40);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a41[32] = {0x8d,0xae,0x9f,0xfb,0x47,0x84,0xf9,0xdd,0xd2,0x68,0x49,0x1f,0x91,0x21,0xa4,0x3e,0xa2,0xb8,0xf9,0xed,0x25,0x80,0x5f,0xa5,0x89,0x35,0xcf,0x07,0x5c,0xc6,0x5e,0x56};
    static std::uint8_t c41[32] = {0x35,0xdb,0xe3,0xe8,0x1d,0x42,0xe7,0x6f,0x9a,0x0e,0x0e,0x0f,0x9c,0xb2,0xb9,0x79,0x7a,0xa4,0xaa,0x1c,0x85,0xc0,0x90,0x66,0x89,0xa7,0x76,0x13,0x61,0x61,0x56,0xf8};
    sm9_fp_from_bytes(a, a41);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a42[32] = {0x62,0x91,0x73,0x47,0x76,0xc4,0xeb,0x43,0xfc,0xf9,0x0d,0xee,0xa9,0x5c,0xfe,0xd6,0xe9,0x35,0xd2,0xcc,0x15,0x73,0x3c,0xd2,0x38,0x05,0x5b,0x51,0xa9,0x0b,0xdd,0x71};
    static std::uint8_t c42[32] = {0x46,0x50,0x23,0xc3,0xf5,0x19,0x40,0x84,0xc6,0x2b,0x4e,0x4f,0xca,0x8d,0x8f,0x3a,0xb6,0x4d,0x2b,0xe8,0xdb,0x0d,0x98,0xbc,0xcd,0x97,0x9e,0x61,0x1d,0xf5,0x60,0x96};
    sm9_fp_from_bytes(a, a42);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a43[32] = {0x60,0xa2,0x21,0x81,0xf6,0xbc,0x93,0x2d,0x01,0x8e,0x68,0x24,0xcd,0xee,0x39,0x2d,0xf0,0x66,0x15,0xf4,0xf7,0xb2,0x21,0xe0,0x0c,0xa7,0x2d,0xab,0xe6,0x12,0x4b,0x52};
    static std::uint8_t c43[32] = {0xaa,0xd9,0x2e,0x23,0xf4,0x39,0x83,0xbf,0xfc,0x5a,0x47,0x02,0xe9,0x01,0x00,0xb4,0x46,0xf7,0x44,0x01,0x9e,0x13,0x1f,0x1b,0xb1,0x1e,0xe6,0x02,0x44,0x6b,0x4a,0x1f};
    sm9_fp_from_bytes(a, a43);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a44[32] = {0x09,0x6c,0xca,0xa8,0x4f,0xfe,0x4f,0xd7,0xa9,0x6f,0x61,0x3b,0x77,0xaa,0x3d,0xe6,0x22,0xca,0x8b,0x97,0x79,0x15,0xae,0x56,0x10,0x66,0x3b,0xcd,0x28,0x2a,0xa8,0x16};
    static std::uint8_t c44[32] = {0x2a,0xfc,0xf5,0x8d,0x46,0x97,0x36,0x27,0x34,0xee,0xca,0x8f,0xb7,0x81,0xa2,0xcb,0x25,0x8e,0x3b,0xfe,0xea,0xfb,0xfe,0x6f,0xfb,0x2a,0xd2,0xea,0xba,0xe3,0xc6,0x05};
    sm9_fp_from_bytes(a, a44);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a45[32] = {0x53,0x1b,0x1e,0xd4,0x71,0x5e,0x73,0x5b,0x37,0x2e,0xa4,0x1e,0x4e,0x10,0xa2,0x3f,0xb0,0x21,0xea,0x9b,0x53,0xb7,0x92,0x5e,0xb0,0x23,0x96,0xb5,0x96,0x5b,0x9b,0xcd};
    static std::uint8_t c45[32] = {0x41,0x9c,0x53,0xdc,0x49,0x12,0x7f,0x91,0x0f,0x73,0x50,0x30,0xd0,0x37,0x5f,0x7f,0x26,0x11,0x1f,0x33,0x10,0x98,0x03,0xbd,0x66,0x74,0x3f,0xce,0xe4,0x3b,0x37,0xbf};
    sm9_fp_from_bytes(a, a45);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a46[32] = {0x57,0x06,0x47,0xeb,0x8d,0xf1,0xbb,0xc0,0xa7,0xac,0xad,0x1c,0x75,0xbc,0x7b,0xaf,0x28,0x0a,0x55,0x16,0xd9,0x3b,0x83,0xec,0x94,0xf8,0x21,0x2a,0x35,0xcb,0xf0,0x95};
    static std::uint8_t c46[32] = {0x0a,0xce,0x6b,0x19,0x8c,0x7d,0x9e,0xd6,0x97,0x16,0x83,0x60,0x39,0xe4,0x2b,0x04,0xca,0x71,0xa3,0x9b,0x73,0xaa,0x6c,0xab,0xac,0x98,0xdb,0x0a,0xf3,0xf3,0xe0,0x4c};
    sm9_fp_from_bytes(a, a46);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a47[32] = {0x06,0x09,0x3f,0xb5,0x47,0x1d,0xb2,0xcc,0x9a,0x64,0xe2,0x3e,0xf3,0xe6,0x49,0x98,0xb3,0x24,0xed,0xa0,0x67,0x09,0xe1,0xe0,0xb1,0xec,0xe4,0x98,0xdb,0x81,0x73,0xaa};
    static std::uint8_t c47[32] = {0xb0,0x10,0x36,0x22,0x06,0xe2,0x02,0x52,0xfd,0x47,0x07,0x6c,0xb5,0xfc,0x7f,0xff,0xce,0x69,0x77,0xee,0xff,0x03,0x51,0x0c,0x14,0xe7,0x29,0x64,0x24,0xe4,0x3c,0x0d};
    sm9_fp_from_bytes(a, a47);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a48[32] = {0x1e,0x88,0xfd,0xef,0x15,0x41,0xab,0xd6,0x7e,0x7a,0x97,0xe6,0x4f,0xcd,0xd4,0xd0,0x38,0x89,0xcb,0x98,0xbb,0x7a,0x4a,0x6c,0xfe,0x1f,0x1d,0x54,0x6e,0x7f,0xf6,0xfd};
    static std::uint8_t c48[32] = {0x1e,0x1e,0xa6,0x89,0xe4,0x10,0xd8,0x79,0x9d,0xfd,0xf1,0xad,0x49,0x36,0xd5,0x91,0xf3,0x04,0x42,0xa7,0x00,0x51,0xee,0xf2,0x56,0xf3,0x56,0x2e,0xe0,0xf8,0xa3,0x6b};
    sm9_fp_from_bytes(a, a48);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a49[32] = {0x5b,0x62,0xaa,0xf1,0x2f,0x1d,0xfb,0x68,0xb3,0x93,0x8f,0x79,0x90,0x19,0x57,0xd0,0x58,0x3f,0xc1,0x1b,0xcf,0x66,0x05,0xb8,0x57,0xde,0x23,0xcb,0x55,0xd2,0xa6,0x36};
    static std::uint8_t c49[32] = {0x76,0x6b,0xd9,0xd8,0x06,0xc4,0xbd,0x92,0xdc,0x6b,0x4a,0xa4,0xdd,0xd3,0xf4,0x8e,0x01,0xff,0x4a,0x46,0x68,0x40,0x7b,0x00,0x75,0x2b,0x72,0xcf,0xb1,0x7a,0xfa,0x96};
    sm9_fp_from_bytes(a, a49);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a50[32] = {0x1a,0x50,0x27,0x61,0x64,0x68,0xe9,0x0e,0xa9,0x92,0x4b,0x8b,0x27,0x5a,0x64,0xac,0x4f,0x52,0x76,0xd3,0xd1,0x62,0x83,0x9d,0x66,0x74,0x28,0x92,0x81,0xac,0xc1,0x3b};
    static std::uint8_t c50[32] = {0x04,0x81,0xbc,0x80,0xa1,0xb9,0xb7,0xf4,0x44,0xa6,0xe9,0xeb,0x17,0x23,0x4d,0x1b,0x06,0x13,0xe2,0x7d,0x88,0xfe,0x9b,0xa0,0xc7,0x39,0x46,0xaa,0x1a,0xc6,0xd8,0xec};
    sm9_fp_from_bytes(a, a50);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a51[32] = {0x3a,0x4e,0xf0,0x33,0xeb,0x0d,0xc5,0x50,0x59,0x15,0x4f,0x7d,0xaa,0x6d,0x8b,0x96,0x2c,0x68,0x4b,0x9b,0x8a,0xb4,0x2f,0x7e,0x92,0x01,0xc2,0xc2,0x39,0x9a,0x8f,0x72};
    static std::uint8_t c51[32] = {0x7d,0xf9,0xd9,0xc7,0x47,0x2a,0x6f,0xf8,0x42,0x6a,0x80,0x48,0x44,0x12,0x95,0xa5,0xfb,0x7c,0x03,0x07,0x9d,0x26,0xe4,0x1e,0xc8,0xb2,0xe7,0xe1,0xaa,0x16,0xcd,0x6e};
    sm9_fp_from_bytes(a, a51);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a52[32] = {0x0b,0xfe,0xe0,0x7d,0x19,0xe0,0x5a,0xba,0xb7,0x14,0xb6,0x5e,0xf2,0x07,0x9e,0x90,0x57,0xf9,0x8f,0x08,0x84,0x6c,0x6e,0x40,0xae,0xa6,0xa6,0xb6,0x3a,0x2f,0xf8,0x92};
    static std::uint8_t c52[32] = {0x44,0xc0,0x27,0x4a,0xe0,0x89,0x6f,0x73,0xec,0x68,0xf1,0x82,0x2f,0x69,0x88,0x40,0x3e,0xb7,0xce,0xbf,0xf8,0x59,0x76,0x55,0x86,0xc8,0x75,0x4a,0x7a,0x2f,0x06,0x92};
    sm9_fp_from_bytes(a, a52);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a53[32] = {0x73,0xf1,0xd3,0xfb,0xc3,0x4d,0x02,0x82,0x16,0x3a,0xf9,0x4e,0x89,0x0f,0x90,0xf7,0xc5,0x13,0xc9,0xe8,0x8e,0xd0,0x39,0x8c,0x51,0x9f,0xce,0xd8,0x66,0x0d,0x84,0x4d};
    static std::uint8_t c53[32] = {0xb4,0xdf,0x32,0x03,0xa5,0x15,0x0c,0xd1,0xd0,0xda,0xfb,0x8d,0x96,0x0e,0x6e,0x5a,0x5c,0x8a,0x79,0xd2,0xe2,0x5d,0xfd,0xa9,0x7b,0xa7,0xc3,0xe5,0x4a,0x46,0xcf,0x6f};
    sm9_fp_from_bytes(a, a53);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a54[32] = {0x5b,0x4e,0xcb,0x41,0x4a,0xa8,0xee,0xc1,0xdc,0x52,0x11,0x93,0xfc,0xa8,0x86,0x1c,0x40,0x09,0x42,0x0e,0x6b,0xbe,0x22,0x9c,0x74,0x98,0x00,0x6b,0x50,0x8f,0xe8,0x72};
    static std::uint8_t c54[32] = {0xa6,0x7d,0x80,0xe0,0x70,0xa8,0xd9,0x01,0x0b,0xc1,0x14,0xfe,0x4e,0x4b,0x71,0xfe,0xe6,0x28,0x9d,0x27,0xca,0x5c,0x66,0x6c,0x2e,0xb8,0x9f,0x7f,0x11,0xff,0xab,0x30};
    sm9_fp_from_bytes(a, a54);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a55[32] = {0x8c,0xbc,0x7a,0xa7,0xec,0x0e,0xf6,0x7f,0x62,0x03,0x76,0x4b,0x6f,0x30,0xec,0x1b,0x9e,0xb0,0x22,0x13,0x77,0xf4,0x2f,0xba,0x53,0x96,0x5d,0xe2,0x14,0x9c,0x7a,0x3f};
    static std::uint8_t c55[32] = {0x39,0x3b,0x3d,0xeb,0xfa,0xa0,0x39,0x14,0xe6,0xd9,0x4d,0xfa,0xd0,0x45,0xdb,0x70,0xf0,0x84,0x16,0xa4,0x17,0xb1,0x55,0x97,0xf5,0x67,0xf3,0x7e,0xa6,0xdd,0xd6,0xe5};
    sm9_fp_from_bytes(a, a55);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a56[32] = {0x7d,0x12,0x3b,0x15,0x08,0x0c,0xa3,0x13,0x1d,0xa6,0x43,0x52,0x38,0x56,0x1c,0xfb,0xf2,0x3e,0xc7,0xc4,0x00,0xc3,0x63,0x07,0xa9,0x42,0x66,0xe1,0x43,0x9b,0x62,0x15};
    static std::uint8_t c56[32] = {0x5e,0x41,0x85,0x11,0x7f,0x3e,0xee,0x5b,0x6c,0xef,0x02,0x04,0x17,0x51,0xcb,0x70,0x3d,0x01,0x90,0xe0,0x4d,0xe0,0x98,0xf4,0xbb,0xab,0xa7,0xed,0x68,0x5c,0x31,0xef};
    sm9_fp_from_bytes(a, a56);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a57[32] = {0x8a,0x13,0x73,0x20,0xf0,0xae,0x45,0xbf,0x40,0xaf,0x36,0x08,0x34,0x4b,0xcc,0x78,0xdd,0xad,0x3d,0xc3,0xdd,0xaf,0x0e,0x03,0x8e,0x98,0xfe,0x71,0xf3,0xef,0xae,0xb9};
    static std::uint8_t c57[32] = {0x19,0x76,0x63,0x2e,0x32,0xfd,0xf3,0x90,0xd4,0x78,0xbd,0xad,0x29,0xe8,0xfb,0xf3,0xfa,0xb8,0x5e,0xa8,0x04,0x6a,0xc3,0x76,0xa0,0x9a,0xcb,0xdb,0x8b,0x64,0xc7,0x00};
    sm9_fp_from_bytes(a, a57);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a58[32] = {0x2d,0x6b,0xf9,0xe2,0x6c,0x09,0x79,0x73,0xc8,0xee,0x56,0x9f,0xc5,0x86,0x5d,0xa1,0x0b,0x1b,0x43,0x0b,0x23,0xed,0xc8,0xbe,0x58,0xbf,0x76,0xe1,0x26,0x57,0x63,0x5c};
    static std::uint8_t c58[32] = {0xac,0xf6,0xf5,0xf8,0x5a,0x6a,0x4d,0x6d,0xae,0x62,0xd2,0x20,0x77,0xf0,0x1d,0x7c,0x35,0xb2,0x60,0xeb,0xb9,0x1f,0xe9,0x76,0x85,0x59,0x7f,0xe4,0x19,0xdb,0x6d,0xde};
    sm9_fp_from_bytes(a, a58);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a59[32] = {0x0a,0xca,0x59,0x2e,0x5b,0xb8,0xb7,0x16,0x12,0x80,0xdf,0xa7,0x82,0xf5,0xc4,0x19,0xd2,0xf0,0x5b,0xaf,0xf2,0x4b,0x9f,0x90,0xec,0x2b,0xee,0xc0,0x25,0x8e,0xe8,0x08};
    static std::uint8_t c59[32] = {0x68,0xa7,0xdd,0xaf,0xc7,0x0e,0x2a,0x35,0x34,0x8c,0xe4,0x45,0xd8,0x6d,0xb0,0x61,0xff,0xa8,0x90,0xf3,0x7b,0x86,0xfe,0x72,0xb6,0xa6,0xa2,0x05,0x2e,0x18,0x49,0x95};
    sm9_fp_from_bytes(a, a59);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a60[32] = {0x30,0xac,0x1e,0x8c,0x8a,0x32,0xc5,0xf8,0x9c,0xe9,0xeb,0x4f,0x0e,0xc1,0x84,0xf4,0x7f,0x5d,0xd3,0x0a,0x55,0x64,0x19,0x9a,0x1f,0xae,0x68,0xdf,0x64,0x32,0x4d,0xd5};
    static std::uint8_t c60[32] = {0x8e,0x38,0x99,0x51,0xd3,0xbe,0x78,0xc7,0xb0,0xc1,0x1f,0x76,0x4a,0x7a,0x34,0x72,0xfe,0x36,0xbc,0xcf,0x7d,0x70,0x5a,0xe2,0xb7,0xb7,0x88,0x24,0x9f,0x7b,0x2d,0xd4};
    sm9_fp_from_bytes(a, a60);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a61[32] = {0x9d,0x8d,0x86,0x39,0xf5,0x6f,0x01,0xca,0xb4,0xd1,0xd0,0x0a,0x0a,0x38,0x0f,0xff,0x4f,0xb1,0x99,0x74,0xfa,0x53,0x93,0x80,0xf6,0xfd,0x5b,0x1d,0xb3,0xbe,0x29,0x59};
    static std::uint8_t c61[32] = {0xaa,0x67,0x4c,0x34,0x15,0x71,0xf4,0x3f,0xf0,0xb0,0x5e,0x27,0x9a,0xaf,0xe3,0x71,0xd8,0x16,0xef,0x82,0x9b,0xd9,0x38,0xd0,0x86,0xa8,0xcf,0x25,0xb2,0x58,0xc0,0x78};
    sm9_fp_from_bytes(a, a61);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a62[32] = {0x78,0xe9,0xa8,0x72,0x15,0x3e,0x81,0x27,0xf1,0xc2,0x5a,0x0b,0xd7,0xb9,0x23,0x7c,0x14,0x76,0x75,0xd6,0x99,0xed,0x0c,0xfb,0xc2,0x97,0xba,0x9c,0x0a,0x4a,0x2a,0x76};
    static std::uint8_t c62[32] = {0x02,0x2d,0x01,0x5e,0x6e,0x30,0xa8,0x2f,0xd4,0xb9,0x58,0x8d,0xdf,0xb1,0x0a,0xbb,0xdf,0x5e,0x8c,0xc7,0x3c,0xf5,0xdb,0x9f,0x7e,0x73,0x82,0x91,0xf7,0xfb,0xa2,0x67};
    sm9_fp_from_bytes(a, a62);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a63[32] = {0x17,0xa4,0x7e,0x42,0x8b,0x75,0x2a,0x2c,0xbf,0xbb,0xe1,0x0e,0x0a,0x5b,0xa8,0xfa,0xab,0xe0,0x50,0xe0,0xd7,0x8c,0x61,0xef,0x65,0x47,0x81,0xc7,0x6a,0xe8,0x2c,0x1d};
    static std::uint8_t c63[32] = {0x10,0x05,0xbd,0x5d,0x48,0x5f,0x4d,0xfc,0x9b,0x66,0xcd,0xde,0x42,0x5e,0xf7,0x4a,0x1a,0x3d,0x81,0x00,0xb3,0xb4,0x7d,0xae,0x23,0xc7,0xe5,0xbc,0xa6,0xae,0x76,0x0b};
    sm9_fp_from_bytes(a, a63);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a64[32] = {0xa8,0x52,0xb8,0x1e,0x86,0x0c,0x7b,0xc3,0x8c,0x36,0x8e,0xc3,0x21,0x2d,0x4e,0xf5,0xb0,0x3e,0xdd,0x30,0xc0,0x1b,0xc6,0x5c,0x28,0x49,0x66,0x88,0xfc,0x93,0x98,0x3c};
    static std::uint8_t c64[32] = {0xb3,0x65,0x2f,0xbe,0xf5,0x1b,0xf6,0xc0,0x3b,0x84,0x0d,0x37,0xd0,0xbb,0xe4,0xd1,0xa9,0x98,0x5b,0xd8,0x47,0xc2,0x3e,0x2e,0xeb,0xae,0xc3,0x9b,0xc5,0xa5,0xa8,0xcf};
    sm9_fp_from_bytes(a, a64);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a65[32] = {0x4b,0xef,0x26,0x33,0x91,0x35,0x04,0x40,0xc5,0x66,0x31,0x89,0x99,0x99,0x8b,0x6c,0x49,0x98,0x1a,0x54,0x9f,0x37,0xcb,0x3e,0xe3,0x12,0x4f,0x8e,0x10,0x8d,0x27,0xcf};
    static std::uint8_t c65[32] = {0xb4,0x23,0xe9,0x52,0x8d,0x0a,0xc4,0x5a,0xa3,0x90,0x03,0xa2,0xa9,0x59,0x69,0xaf,0x2f,0x57,0x8b,0x1c,0x35,0xa6,0x10,0x1a,0xba,0xb2,0xcc,0xb1,0xc5,0xbe,0x43,0xcb};
    sm9_fp_from_bytes(a, a65);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a66[32] = {0x0d,0x8a,0x82,0xae,0xd8,0xda,0x88,0x64,0xf0,0xdd,0xc9,0xb4,0xbc,0xfd,0xe7,0xe4,0x1d,0x0a,0x64,0x51,0x54,0x25,0xe2,0x9c,0xff,0x18,0xec,0x51,0xb2,0x9c,0x51,0xa6};
    static std::uint8_t c66[32] = {0xad,0xd4,0xe0,0x73,0xff,0x3b,0x5d,0x98,0xa1,0x0c,0xe6,0xe6,0x22,0x4e,0x94,0x50,0x3e,0xde,0x35,0x6a,0x7f,0x4f,0xea,0xa7,0xc6,0x66,0x36,0xd1,0x0e,0x1d,0x87,0xb9};
    sm9_fp_from_bytes(a, a66);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a67[32] = {0x4d,0xe9,0xb6,0xb4,0xd1,0xe2,0x81,0xd7,0xc0,0x39,0xb1,0x15,0x67,0x02,0x17,0xc1,0x67,0x18,0x02,0xc3,0x7c,0x7e,0x87,0xdd,0xfb,0xe1,0x8c,0x7c,0x47,0x22,0xdd,0x63};
    static std::uint8_t c67[32] = {0x60,0x8d,0x15,0x3c,0x7c,0x17,0xc2,0x15,0xd8,0xc9,0x81,0x26,0x4d,0xd2,0x8b,0x55,0x11,0x43,0x52,0x9e,0x82,0x78,0x55,0x3d,0x22,0x1b,0x9e,0xf5,0x3b,0x34,0x90,0x4f};
    sm9_fp_from_bytes(a, a67);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a68[32] = {0x69,0xf9,0x16,0x6c,0x36,0xad,0x64,0x21,0xa4,0xb9,0x3e,0xa5,0x6f,0x29,0x41,0xee,0xa0,0xce,0xa4,0x45,0x8d,0x2b,0x51,0xef,0x45,0x93,0x40,0x24,0xfa,0xee,0x65,0x93};
    static std::uint8_t c68[32] = {0x6d,0xa8,0x24,0x63,0xf2,0x97,0x9e,0xcb,0x91,0xf0,0x33,0x84,0x60,0x9c,0x61,0x73,0x12,0x8b,0xe8,0x62,0x88,0x83,0x13,0xd9,0x17,0xe5,0x24,0xbd,0x35,0x0d,0xb8,0x84};
    sm9_fp_from_bytes(a, a68);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a69[32] = {0x14,0x67,0xa1,0x16,0xcf,0x9f,0x5e,0x1d,0x84,0x8f,0x51,0x2b,0xc1,0xc8,0xb6,0x3d,0x06,0xd3,0x41,0x01,0x95,0xef,0x65,0x6e,0xf7,0xcf,0x9e,0x31,0x84,0xa1,0x58,0x82};
    static std::uint8_t c69[32] = {0x88,0x5c,0x28,0x04,0x30,0xab,0x5e,0x85,0x88,0x77,0xa5,0x4c,0x70,0x0b,0x63,0xae,0x9a,0xd0,0x10,0xa8,0x62,0xcb,0xd0,0x05,0x24,0xdf,0xce,0xd8,0x17,0x2d,0x33,0xc1};
    sm9_fp_from_bytes(a, a69);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a70[32] = {0x89,0xa3,0xf5,0xcd,0xa9,0x25,0x84,0xdb,0xea,0xc9,0x35,0x66,0xc9,0xcf,0x9d,0x59,0x7c,0xa5,0x07,0x59,0xb3,0x5e,0x6a,0x0a,0x34,0x36,0x6c,0x91,0x54,0x0d,0x1e,0x8e};
    static std::uint8_t c70[32] = {0x26,0x49,0x3c,0x6a,0x14,0xa1,0x91,0x58,0xad,0xa6,0xa4,0x19,0xc0,0xd4,0x20,0xff,0x94,0xa5,0x69,0xb1,0x92,0x97,0x29,0x67,0x5c,0xf6,0xaa,0x57,0x66,0x3a,0x20,0xd6};
    sm9_fp_from_bytes(a, a70);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a71[32] = {0x4c,0x10,0x2a,0x68,0x3d,0x79,0x9d,0xcf,0xe4,0x0b,0x63,0x5a,0x45,0xff,0x9a,0x45,0x78,0x94,0xfb,0xd2,0x42,0x68,0x1a,0x28,0x98,0x27,0xfe,0x38,0x30,0x2b,0x13,0x2e};
    static std::uint8_t c71[32] = {0x52,0xfc,0x8f,0x1b,0x11,0x1d,0x88,0x31,0xce,0x30,0x56,0xf6,0x06,0xf2,0x24,0xc2,0xc1,0x80,0xe4,0x60,0x7e,0xf5,0xd1,0x8c,0xb8,0xad,0x25,0xb1,0x20,0xc6,0xbe,0xd3};
    sm9_fp_from_bytes(a, a71);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a72[32] = {0x81,0x45,0xd8,0x20,0xdc,0xf1,0xc2,0x6a,0x9a,0x8e,0xd7,0x34,0x84,0x41,0x30,0xf2,0xb7,0xb5,0x48,0x80,0xb3,0xe8,0xf5,0x2b,0x06,0x3d,0x1e,0xf6,0x1d,0xdf,0x99,0x1f};
    static std::uint8_t c72[32] = {0x54,0x6c,0xfc,0xd3,0x2d,0xb1,0xab,0x10,0x75,0xa6,0xda,0x9e,0xb6,0x42,0x8f,0x66,0x66,0x1e,0x71,0x64,0x65,0x4d,0x9c,0xd2,0x90,0x1c,0x67,0x51,0x2b,0xa1,0xfc,0x97};
    sm9_fp_from_bytes(a, a72);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a73[32] = {0x38,0xe3,0x79,0x75,0xa8,0xa2,0xd9,0x45,0x01,0x53,0xd0,0xfd,0x9d,0xa5,0xcc,0xcc,0x7e,0xc5,0x26,0x08,0xfd,0x49,0xec,0x90,0xe7,0x81,0x9b,0xe8,0xe4,0x80,0x0c,0x94};
    static std::uint8_t c73[32] = {0x36,0xfd,0xf0,0xde,0x86,0xd9,0xd3,0xdc,0x68,0x64,0xd1,0x8c,0x22,0x1a,0xd3,0x23,0xde,0xa9,0x88,0xab,0xbd,0xb5,0xf7,0x06,0x61,0xf8,0x41,0x0f,0x47,0xed,0x66,0xf8};
    sm9_fp_from_bytes(a, a73);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a74[32] = {0xa7,0x75,0x3e,0xdc,0xe3,0x9e,0x28,0xc3,0xf8,0x01,0x37,0x27,0x3b,0xda,0x8f,0x95,0x57,0x6a,0x02,0x6e,0xdb,0xab,0x50,0x5e,0xbf,0xf4,0x40,0x49,0xa5,0x51,0x2c,0x39};
    static std::uint8_t c74[32] = {0xb2,0xd0,0x7e,0xbc,0x92,0xb6,0x3a,0x8a,0x29,0xc3,0x88,0x56,0x04,0x39,0xe0,0x6e,0xfb,0xc3,0xbf,0x9c,0x86,0x28,0xbd,0x55,0xcc,0xa7,0x97,0x48,0x56,0x79,0x5d,0x2f};
    sm9_fp_from_bytes(a, a74);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a75[32] = {0x57,0x1c,0x09,0x1e,0xa2,0xde,0x70,0x57,0x9c,0x49,0x4a,0xcb,0xb7,0xf4,0x51,0xfe,0x27,0x6e,0xa9,0x97,0xa9,0x0c,0x05,0xc2,0x9f,0x4a,0x4f,0xbc,0x8b,0x58,0xcd,0xb7};
    static std::uint8_t c75[32] = {0x82,0x1f,0x18,0x64,0xa4,0x05,0x5b,0x23,0xa1,0xab,0x0e,0x9f,0x0d,0x10,0xf9,0x80,0x43,0x46,0x03,0xca,0xed,0x11,0x64,0xbb,0xf3,0xe3,0x14,0x72,0xd8,0x89,0xd8,0x4c};
    sm9_fp_from_bytes(a, a75);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a76[32] = {0x00,0xcb,0xa8,0x16,0x86,0x5c,0xc8,0xcf,0x83,0x1e,0xc5,0x52,0xc7,0x7a,0x6c,0xc2,0x1b,0xaf,0x3e,0xeb,0xa6,0xdf,0x30,0xf4,0xf2,0xea,0x24,0xb3,0x5c,0xa1,0xaa,0xe8};
    static std::uint8_t c76[32] = {0x91,0x9b,0xf0,0x87,0xaa,0x2a,0x9b,0x7a,0x33,0x90,0x4c,0xe2,0xa4,0x4e,0xd5,0x46,0x22,0xe4,0x05,0xdd,0x74,0x4e,0x8a,0x60,0xe4,0xde,0x25,0x90,0x3a,0x2f,0x5d,0x3d};
    sm9_fp_from_bytes(a, a76);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a77[32] = {0x7b,0xff,0xbe,0x9c,0x84,0x7c,0x16,0x94,0x6b,0x15,0xd4,0xe7,0x2d,0x82,0x4b,0x10,0x86,0x31,0xb1,0x90,0xed,0x5d,0xc9,0x1e,0x63,0x2f,0x11,0xad,0xa8,0x3a,0xd2,0x97};
    static std::uint8_t c77[32] = {0x3c,0xe1,0xe2,0x8e,0x8f,0xf1,0x8e,0x00,0xcf,0x3f,0x54,0xef,0x69,0xde,0x04,0x80,0xd1,0x2a,0xff,0xe6,0x4a,0xc1,0xb8,0x4e,0xba,0xfd,0x1e,0x61,0x48,0x79,0x9e,0x70};
    sm9_fp_from_bytes(a, a77);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a78[32] = {0x92,0x75,0x3c,0x98,0xfb,0x4f,0xc9,0xe4,0xf8,0xa2,0xb2,0x2d,0x48,0x76,0x70,0x7e,0x62,0x92,0x9e,0xe4,0x11,0x70,0xbe,0x9a,0xf8,0x47,0x54,0xe0,0xd1,0xa2,0xbf,0xdf};
    static std::uint8_t c78[32] = {0x67,0xe8,0xde,0xd8,0x15,0x83,0xa2,0x8c,0x2b,0x5f,0x0a,0x2e,0x65,0xfb,0x1c,0x7a,0x68,0xda,0x96,0x0e,0xea,0x72,0xd9,0x40,0x65,0xdb,0x11,0x84,0xf6,0x4e,0xa5,0xca};
    sm9_fp_from_bytes(a, a78);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a79[32] = {0x1d,0x13,0x99,0xfc,0xe8,0x17,0x25,0x69,0xa2,0xb4,0xae,0x06,0x39,0xd2,0x0f,0xee,0xce,0x43,0x29,0xb4,0xa0,0x38,0x7d,0x0a,0xf4,0x29,0x1a,0xc9,0x10,0x49,0xa7,0xfa};
    static std::uint8_t c79[32] = {0x03,0xf2,0xea,0xd7,0xf7,0xf4,0x64,0xfe,0x14,0xe7,0xf1,0x87,0x10,0xcd,0x97,0xdf,0x2b,0xe0,0xdd,0x06,0xde,0xc0,0xd9,0x5d,0xdd,0x30,0x39,0xaf,0x04,0x95,0xeb,0xdc};
    sm9_fp_from_bytes(a, a79);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a80[32] = {0x3c,0x2f,0x0e,0x76,0x70,0x1b,0x88,0x94,0xfe,0xb7,0x76,0x4a,0xed,0xa6,0xf4,0xd9,0x7f,0xec,0x38,0x81,0x94,0x1d,0x98,0x31,0x34,0x45,0x0a,0xd8,0x78,0x20,0x54,0xdd};
    static std::uint8_t c80[32] = {0x95,0xb1,0x0c,0xc7,0xb3,0x46,0x44,0x84,0x94,0x13,0x24,0x64,0xda,0x71,0x7d,0xcd,0x5e,0x6d,0xe9,0xab,0x8d,0x22,0x12,0x11,0x3d,0x72,0xa8,0xb5,0x8e,0xbd,0x78,0x67};
    sm9_fp_from_bytes(a, a80);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a81[32] = {0x3c,0x01,0xbc,0x2b,0xb2,0x49,0xda,0xb5,0x87,0x28,0xe5,0x28,0xad,0x54,0x7c,0xd8,0x09,0xa2,0xbf,0x2e,0x3b,0x24,0xcd,0xc5,0x0a,0x5c,0x1e,0x50,0x80,0xfb,0x30,0xf3};
    static std::uint8_t c81[32] = {0x0e,0xbb,0xe6,0x26,0xf6,0x5e,0x20,0xfc,0xf5,0x5e,0x96,0xb8,0x26,0x64,0xc3,0xf7,0x35,0xf9,0x9e,0xe9,0x4f,0x6f,0xb8,0x20,0xf5,0x7e,0x4b,0xd9,0xb5,0x24,0x20,0x00};
    sm9_fp_from_bytes(a, a81);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a82[32] = {0x9c,0x46,0xec,0x75,0xe9,0x44,0x50,0xfa,0x1b,0x99,0x53,0xff,0xf3,0x71,0xf4,0x01,0x59,0x45,0x50,0xdf,0x07,0x6a,0xf4,0x19,0x96,0xf9,0x9c,0x34,0xab,0xd4,0x52,0x89};
    static std::uint8_t c82[32] = {0x44,0x96,0xcc,0x18,0x91,0x14,0xa9,0x82,0x33,0x88,0xeb,0xb8,0xbc,0x6d,0x25,0x8b,0x32,0x77,0xfa,0x33,0xf7,0x5a,0xc6,0xe7,0x36,0x0a,0x83,0x31,0x7a,0x12,0x36,0xa8};
    sm9_fp_from_bytes(a, a82);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a83[32] = {0xa5,0x53,0xfc,0x7a,0xb5,0x59,0x98,0x30,0x24,0x02,0x9d,0xdf,0x3a,0xb7,0x84,0x95,0xd2,0x9c,0x7e,0x1f,0xf6,0x88,0xf2,0x8a,0x29,0x98,0xd2,0x07,0xbb,0xaa,0xaa,0x82};
    static std::uint8_t c83[32] = {0x3a,0xbe,0x21,0x3d,0xce,0x9e,0x3e,0xfa,0x87,0xbc,0x25,0x7e,0x45,0x96,0xcd,0x13,0x0d,0x89,0x5b,0x1b,0x07,0xe6,0x0f,0x11,0xbb,0xa2,0x21,0xcb,0xdc,0x23,0x7f,0xce};
    sm9_fp_from_bytes(a, a83);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a84[32] = {0x0b,0x7c,0x38,0xd9,0x15,0xc5,0x36,0xcf,0x7e,0x12,0x75,0x9a,0x8b,0xdf,0x8a,0x58,0xa8,0x4a,0x97,0x81,0xc1,0x02,0x9f,0xb3,0x80,0x7d,0xf4,0xd8,0x7e,0xe5,0xde,0xef};
    static std::uint8_t c84[32] = {0x01,0xc7,0x09,0x1b,0xe1,0xc4,0x27,0xb1,0xec,0x0d,0x72,0x3b,0x81,0xf6,0x28,0x4d,0xdb,0xcb,0xe7,0xcd,0xce,0xff,0x93,0xb0,0x54,0xc1,0x0e,0xa9,0x14,0x71,0x7f,0x07};
    sm9_fp_from_bytes(a, a84);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a85[32] = {0x82,0xa9,0xc0,0x08,0x1e,0xd8,0xdd,0x89,0xdc,0x8e,0x86,0x0f,0x46,0x75,0x8d,0x61,0xc4,0x7c,0x7d,0xa2,0x09,0x07,0x71,0x0a,0x44,0x4a,0x7c,0x02,0xfc,0x5c,0x89,0xe3};
    static std::uint8_t c85[32] = {0x83,0x75,0x96,0x72,0x26,0x36,0x79,0x67,0xef,0x43,0xc1,0x4f,0xf1,0x53,0x08,0xd4,0x68,0xca,0xf3,0x3d,0x1e,0x88,0xe9,0xc7,0x20,0xb9,0xc1,0x92,0xea,0x78,0xfa,0xaf};
    sm9_fp_from_bytes(a, a85);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a86[32] = {0x65,0x06,0xae,0xee,0x0b,0x56,0x2d,0x13,0xcb,0x57,0xd2,0xad,0x57,0x97,0x6f,0xf7,0x1c,0xfb,0x4c,0x1d,0x24,0x0c,0x2f,0x47,0x08,0x23,0xa4,0xf1,0x5d,0xe0,0xcb,0xd2};
    static std::uint8_t c86[32] = {0x54,0x8e,0xd0,0x80,0x35,0x15,0x5c,0x4d,0x2e,0x13,0xd0,0x59,0x7f,0xca,0x30,0x0d,0xbc,0x51,0x11,0x43,0x9f,0x97,0x0f,0x98,0x65,0xa1,0x44,0x0e,0x52,0xac,0xf4,0x20};
    sm9_fp_from_bytes(a, a86);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a87[32] = {0x75,0xf8,0x89,0xe2,0x2e,0x94,0xe2,0x08,0xf7,0x24,0x83,0x57,0xde,0x6f,0x0d,0xd8,0x07,0x1e,0x27,0x6b,0xa8,0x5d,0x83,0x89,0xfe,0xd2,0x75,0x8f,0xca,0x07,0x1c,0xe4};
    static std::uint8_t c87[32] = {0x94,0xbd,0x24,0xe6,0x19,0x06,0x11,0x76,0x9b,0x2b,0x05,0xe5,0xc6,0x6e,0x99,0x06,0x2c,0x3e,0xeb,0x4e,0xee,0x7b,0xa4,0x98,0x3c,0x80,0x72,0x07,0xfc,0x6f,0x22,0x45};
    sm9_fp_from_bytes(a, a87);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a88[32] = {0x06,0x8c,0x58,0xb6,0x0e,0x6c,0x2b,0x80,0x56,0x26,0x1f,0x12,0x97,0xa3,0xe1,0x0f,0xc8,0x3c,0x63,0x0c,0x46,0x53,0x10,0xd2,0x6b,0xfe,0x10,0xda,0x0e,0xcf,0xda,0x7d};
    static std::uint8_t c88[32] = {0x30,0x2a,0x35,0x27,0xca,0x00,0x77,0xf4,0xbe,0x78,0x8b,0xb4,0x35,0x58,0xf7,0xf2,0x38,0x8e,0xcb,0xe2,0x51,0x15,0xb3,0xff,0x62,0x0c,0x88,0xf3,0x89,0x8d,0x45,0x8e};
    sm9_fp_from_bytes(a, a88);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a89[32] = {0x84,0x60,0x97,0xe9,0x8b,0x88,0x85,0x06,0xfd,0x30,0x9e,0xdc,0x38,0x8f,0xb8,0x98,0x06,0x80,0x32,0x56,0x2a,0x7d,0x1e,0x2a,0xa6,0x83,0x57,0xaa,0x01,0x48,0x9c,0xd3};
    static std::uint8_t c89[32] = {0xa1,0xf9,0xf3,0xa2,0x74,0x8a,0x86,0x83,0x9d,0x20,0x44,0x87,0x3f,0x06,0x39,0xb9,0x25,0x99,0xdf,0xa1,0x86,0xf7,0xcf,0x4e,0x2a,0x38,0xe2,0x4c,0x1b,0x65,0x82,0xe1};
    sm9_fp_from_bytes(a, a89);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a90[32] = {0x6c,0xb4,0xed,0x1b,0x57,0x2a,0xd2,0xbc,0x91,0x3b,0xf5,0x49,0x14,0x04,0x78,0x3e,0xe9,0x66,0xf3,0x51,0xb0,0xa1,0x60,0x52,0xce,0x99,0x95,0xd5,0xba,0x31,0x79,0xdf};
    static std::uint8_t c90[32] = {0x5f,0xdc,0x2c,0x38,0x8e,0x71,0x7a,0x13,0xde,0x83,0x15,0x8d,0x0a,0xbe,0x92,0x45,0x6b,0x37,0x9b,0xe1,0xfc,0xe2,0xbf,0xf3,0x91,0xc9,0xf1,0x25,0x10,0x39,0x64,0x0e};
    sm9_fp_from_bytes(a, a90);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a91[32] = {0x72,0x1b,0x66,0xbb,0x0b,0x1f,0xd0,0xaa,0xba,0x83,0x45,0xf5,0x6f,0xa9,0xea,0x57,0x0a,0xcb,0x52,0x91,0x17,0xf1,0x03,0xda,0x99,0x2d,0x18,0x19,0x84,0xbe,0xcb,0xac};
    static std::uint8_t c91[32] = {0x0d,0xdf,0xf0,0x88,0xea,0x54,0xf1,0xc4,0xdf,0x0d,0xf5,0xb0,0xfe,0x01,0xd6,0x1d,0x2f,0x38,0x5a,0xd5,0x13,0x26,0x62,0xd2,0xf3,0x71,0xdc,0x6c,0xe1,0x8e,0xd8,0x91};
    sm9_fp_from_bytes(a, a91);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a92[32] = {0x03,0xc3,0xfe,0x9b,0x1e,0x43,0x62,0x4e,0xbe,0x75,0x22,0x7e,0x3a,0x9f,0xdd,0x3c,0xc8,0x4d,0xfb,0x46,0x26,0x4b,0x7f,0xf6,0x92,0xaf,0xfe,0xe5,0x93,0xc4,0x0b,0x96};
    static std::uint8_t c92[32] = {0x06,0x9f,0xa7,0xcf,0x9c,0x3e,0x9f,0x97,0xdc,0xd9,0x4f,0xf0,0xa2,0x16,0xdd,0x84,0x06,0xba,0x4f,0x8f,0x37,0x69,0xeb,0x35,0x48,0x82,0x61,0x72,0x9f,0xe1,0x18,0x15};
    sm9_fp_from_bytes(a, a92);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a93[32] = {0x85,0x11,0xd7,0xd9,0x93,0x3f,0x62,0x7c,0x77,0x9e,0x9a,0xa3,0x2d,0x74,0x7e,0xcd,0x3c,0x04,0x22,0x8a,0xd1,0x2c,0x12,0x2f,0x01,0x9b,0x23,0xa0,0xd0,0xbc,0x17,0xde};
    static std::uint8_t c93[32] = {0x6c,0x22,0x61,0x0c,0x9e,0x09,0xef,0xd3,0x6e,0x7a,0x34,0x33,0x16,0xcf,0x55,0x1d,0x62,0xc2,0x45,0xc1,0x46,0xdf,0xd0,0x8b,0xbe,0xc4,0x1b,0xe0,0x78,0x53,0x50,0x5b};
    sm9_fp_from_bytes(a, a93);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a94[32] = {0x3b,0xf2,0xf9,0xa3,0xcd,0x0b,0x13,0xe8,0x07,0x3b,0xa6,0xcf,0x58,0x88,0xf5,0xaa,0x6b,0x86,0xc6,0x99,0xd1,0xde,0x2a,0x54,0x4d,0x0b,0x7e,0x75,0xb0,0x9b,0xf7,0x82};
    static std::uint8_t c94[32] = {0x03,0x49,0x2a,0xf1,0x34,0xac,0x21,0xd0,0xac,0x2b,0xab,0x8a,0xbf,0x6d,0x0d,0x14,0x2f,0xa9,0x4a,0xa0,0xac,0x61,0xe4,0x2e,0xab,0x65,0x5b,0x03,0xb4,0x2b,0x4a,0xca};
    sm9_fp_from_bytes(a, a94);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a95[32] = {0x47,0x06,0xec,0x7f,0x5e,0x58,0xdf,0xb5,0x92,0xa2,0xd2,0x4b,0xd0,0xa7,0xee,0xf3,0xd5,0x09,0x07,0xd4,0x3c,0x3d,0x23,0xcc,0xb5,0x38,0xb9,0xb4,0x52,0x74,0x9d,0xe7};
    static std::uint8_t c95[32] = {0x26,0x24,0x65,0xf9,0x99,0x58,0x2c,0x62,0xd8,0x00,0x21,0x1a,0x91,0xd5,0xab,0x47,0xa4,0x88,0x12,0xac,0x80,0x99,0x38,0x00,0x3d,0xbe,0x7b,0xef,0xa5,0xe0,0x3e,0x04};
    sm9_fp_from_bytes(a, a95);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a96[32] = {0x96,0x44,0x35,0x2b,0x92,0x3d,0xa0,0xb3,0xe3,0xd5,0xa9,0xac,0x29,0x4c,0xb5,0x98,0xeb,0x9e,0xf2,0xc9,0xc3,0x99,0x0e,0xa1,0xb2,0xc9,0xb4,0x9f,0xf1,0x10,0xb8,0x89};
    static std::uint8_t c96[32] = {0x1e,0x72,0x61,0xe7,0xb3,0xf7,0x72,0x73,0xa5,0xfa,0x8c,0x5c,0xd7,0xe5,0x47,0x10,0x8c,0x61,0x3e,0xb0,0x41,0xfb,0x5c,0x95,0xf3,0xac,0x49,0xf4,0x1b,0x9a,0x6d,0x96};
    sm9_fp_from_bytes(a, a96);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a97[32] = {0x18,0x11,0x25,0xbd,0xb6,0xe3,0x16,0x32,0x09,0x07,0x17,0x67,0x3d,0x7c,0xb4,0xab,0x1d,0x6a,0x7c,0xdb,0x4b,0x4b,0xd9,0x9e,0x06,0xb3,0xad,0xf8,0x8c,0x68,0x94,0x75};
    static std::uint8_t c97[32] = {0x12,0xfb,0x0e,0x2f,0x98,0x65,0xdf,0x03,0xbb,0x1a,0xfa,0x34,0xec,0x57,0x16,0xbf,0xd7,0x5a,0xc5,0x4b,0x3c,0xb7,0x2c,0x1f,0x49,0x26,0x8b,0x4a,0xe3,0x65,0xf4,0x9a};
    sm9_fp_from_bytes(a, a97);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a98[32] = {0x30,0x46,0xe7,0xcf,0x35,0x7b,0x86,0x05,0xdd,0x9a,0xce,0xdc,0xe3,0xa3,0x7a,0x11,0x16,0xb0,0xa4,0x6e,0x57,0x77,0xb4,0x51,0x93,0x3b,0x03,0xb1,0xf3,0x77,0x7e,0x33};
    static std::uint8_t c98[32] = {0xb1,0xd5,0xb9,0x16,0x5c,0x9d,0xec,0x32,0x45,0x26,0x31,0xc5,0x65,0x01,0x2c,0x48,0xe0,0x42,0x9d,0x2b,0x62,0x9a,0x60,0x79,0xbf,0x7a,0x62,0xcd,0x5b,0xa8,0x54,0xd9};
    sm9_fp_from_bytes(a, a98);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
    static std::uint8_t a99[32] = {0x90,0x7f,0x96,0xea,0x8e,0xbc,0xa7,0x0a,0x58,0xf2,0x99,0xd7,0x7d,0xb5,0x81,0x88,0xb9,0xc0,0x5e,0xac,0x18,0x64,0xd1,0x2e,0x90,0xa9,0x36,0xd0,0xd0,0x2f,0x0d,0x0b};
    static std::uint8_t c99[32] = {0xa1,0x10,0x6f,0x02,0x35,0xb3,0xa1,0x1f,0xd4,0xa7,0x02,0x21,0x48,0xc0,0x8c,0xd3,0x2d,0x02,0xf0,0xf3,0x86,0x89,0x9d,0xc3,0x51,0x7f,0xfa,0xac,0xa9,0x10,0xdb,0x68};
    sm9_fp_from_bytes(a, a99);
    if(sm9_fp_sqrt(c, a))
    {
        throw std::runtime_error("err in call sm9_fp_sqrt");
    }
    sm9_fp_sqr(c, c);
    if (!sm9_fp_equal(c, a))
    {
        throw std::runtime_error("err in sm9_fp_sqrt");
    }
}